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Abstract 



The fields generated by "quarks" in the adjoint representation of SU(2) color 
are analyzed in the scaling region of the four- dimensional lattice theory. Evi- 
dence of vacuum screening of adjoint quarks is obtained from a comparison of 
quark-antiquark (QQ) flux-tubes for quarks in the adjoint ("isospin" j = 1) 
and fundamental (j = 1/2) representations. The component Sj of the color- 
electric field strength in the direction parallel to the QQ axis is calculated. 
Near the quarks the ratio of fields £?=i /£?=i/2 approaches the value 8/3, which 
is equal to the ratio of SU(2) Casimirs. In between the quarks, the ratio falls 
well below 8/3 at large R. Sj also falls off rapidly as a function of distance xj_ 
perpendicular to the QQ axis. However, the ratio Sj = \/ Sj = i/ 2 depends very 
weakly on xj_. The flux-tubes in the two representations thus appear to have 
very similar cross-sections. This result could imply that the QCD vacuum is 
dual to a type I superconductor. 
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Introduction Linear confinement in QCD has been well established by lattice simula- 
tions of the quark-antiquark (QQ) potential [f]. However the physical mechanism underlying 
confinement has yet to be established, and this continues to represent a fundamental chal- 
lenge to our understanding of QCD. 

Valuable information has been obtained from studies of the color fields generated by a QQ 
pair, as this provides a more detailed probe of the system than is obtained from the potential. 
For example, lattice simulations have demonstrated flux-tube formation [2-5]: color fields 
generated by quarks are found to be "squeezed" into a cylinder whose cross-section is roughly 
independent of their separation R. Flux-tube formation provides an attractive description 
of confinement in QCD, and is the basis of many phenomenological models [6]. 

This paper presents the first calculation of flux-tubes generated by adjoint "quarks" in 
the scaling region of four- dimensional SU(2) lattice gauge theory. It is widely expected 
that adjoint quarks should exhibit dynamics with some important differences from quarks 
in the fundamental representation, thus providing a unique probe of confinement physics 
[7-10,4,11]. In fact, there should be no confinement of adjoint quarks: at large separations 
R it should be energetically favorable for the flux-tube to fission, with the formation of a 
pair of color-neutral quark-gluon bound states ("gluelumps") [8,9]. 

Recently, Michael has obtained some indication that the potential for adjoint quarks may 
saturate at very large i?, where the energy in the flux-tube exceeds twice the energy of the 
gluelump [9] (see also Ref. [11]). Further evidence of adjoint quark screening is obtained 
here by comparing the distribution of fields generated by quarks Qj in both the adjoint 
("isospin" j = 1) and fundamental (j = 1/2) representations of SU(2) color. A number of 
highly efficient techniques are exploited to maximize the quality of the results. 

Method The SU(2) lattice theory with standard Wilson action is used. An efficient 
scheme for the introduction of a static Q 3 Q rj pair is to use a nonlocal (or "fuzzy") Wilson 
loop, one that is "spread out" over an extended region that is, in principle, related to the 
physical dimensions of the state [12,13,8,9,1,5]. Given a set of links U^x), a corresponding 
set of "fuzzy" links U*(x) is constructed according to [13]: 

cU,{x)+ UAx)U,(x + v)Ul(x + ji) , (1) 

where c is a positive constant, and Af is an arbitrary normalization (conveniently chosen 
so that det U 1 = 1). This procedure can be iterated. The number of iterations and the 
parameter c are chosen "empirically" so as to maximize the overlap of a fuzzy operator with 
the state of interest. In order to preserve a transfer matrix fuzzing is applied only to spatial 
links (// ^ 4), which are mixed only with their spatial "staples" [13]. Gauge-invariant (fuzzy) 
Wilson loops are built in the usual way: 

Wi(R,T) = ^yTr jn W] j , (2) 

where L denotes the closed loop, and D^Ui] is an appropriate irreducible representation of 
the link ("fuzzy" spatial links and "ordinary" time-like links). W\ can be expressed in terms 
of W 1/2 , W t = (AW? /2 - l)/3. 



U^ 4 (x)=Af 
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Lattice measurements of the QjQj field strengths are obtained from plaquette correlators. 
Results are reported here for the component of the color-electric field in the direction x 
parallel to the QjQj axis 



£ 3 {x) 



P 



141 



14 



(3) 



where U\± is a plaquette with sides along the x and time-like axes. In the continuum limit 
£j reduces to the expectation value of \ Y^ a (E a ■ x) 2 in the presence of a Q 3 Q 3 pair, after 
vacuum subtraction [2]. In most phenomenological models the structure of the flux-tube is 
determined mainly by E a ■ x. This was shown explicitly for j = 1/2, 1 and 3/2 flux-tubes in 
three-dimensional SU(2) [4]. Lattice calculations also indicate that the potential for j = 1/2 
quarks in four dimensions comes mainly from this field [2] . 

An analytical integration is made on links which occur linearly in the observable (in this 



case, on the time-like links in W 3 j 



\dUx\Vi\Ux\t 



[14]: 
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where S is the action, and k{V\ is equal to the sum of the six "staples" coupling to the link 
U\ (det Vi = 1 and k\ > 0). An analytical integration is also performed over the four links 
in a plaquette [15]: 
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where 
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and Z'{^) = dZ(~f)/d~{. The matrix jiJi is equal to the sum of the five "staples" connected 
to the link Ui, excluding the staple which forms part of the plaquette U\u\U2U\ (det J; = 1, 
and ji > 0). Xn(X) is the character of the SU(2) matrix X in the representation of dimension 
n. The characters are generated from the recursion relation Xn = Xn-iX2 ~ X™-2, starting 
from xi = 1 an d X2(X) = ^TrX [2]. The series in Eq. (6) converges very rapidly (at f3 = 2.4 
less than ten terms suffice in practice). An efficient downwards recursion technique is used 
to generate the Bessel functions I n . 

Equation (5) provides for a reduction in errors in Ej by a factor of as much five. It is also 
advantageous to use (Wj) (U) ~ (WjU(xR)) in Eq. (3) [2], where the reference point xr was 
taken to be 8 lattice units in each spatial direction relative to the center of the Wilson loop 
(which for present purposes is sufficient to ensure a good approximation, within statistical 
errors). This tends to reduce the effect of fluctuations in the Wilson loop [2]. 

Results The calculations were done on a 16 4 lattice at f3 = 2.4. More than 10,000 
sweeps were used for thermalization. 1,800 measurements were made of fuzzy Wilson loops 
of sizes RxT from 3 X 1 to 8 X 5. 6,000 further measurements were then made of loops with 
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R = 7 and 8. 100 heat bath sweeps were made between measurements, yielding integrated 
autocorrelation times T int ~ 0.5. Estimates of the statistical errors were obtained using 
the jackknife method. However, measurements of different observables tend to be strongly 
correlated if they were measured simultaneously on a given lattice. 

Twenty iterations of the fuzzing procedure Eq. (1) were used, with c = 2.5. Figure 1 
shows some results for the "time-dependent" estimates of the potential: 



V estj (i?,T) = - In 



Wi(R,T-l) 



(7) 



V es t ,j(R, T) becomes independent of T for T > 3, within statistical errors. The same is true 
of £j(R, T), although the statistical errors in the fields are much larger than in the potential 
(useful results for R = 7 and 8 could only be obtained for T = 3). Figure 2 shows the 
time-dependence of the fields at the center of the Wilson loop for R = 3 and 5. 

Figure 3 shows evidence of screening in the adjoint potential that becomes more pro- 
nounced as R is increased. At R x T = 8 x 4, for example, Kst,i/Kst,i/2 = 2.47 ± 0.02. 
Michael has obtained results for larger R [9]; at R = 12, V1/V1/2 is between about 2.1 and 
2.3 (the uncertainty comes from an estimate of the systematic errors in the extrapolation 
to T — > 00). These results are to be compared with the ratio of Casimirs j(j + 1) of the 
two representations, equal to 8/3 = 2.66 (the Casimir can be interpreted as the squared 
"charge" of the quark [6]). 

Screening more evident in the field strength. Figures 4 and 5 show the fields as functions 
of position x\\ in the plane of the Wilson loop, using data at T = 3. The quarks are located 
at X|| = ±i?/2. For odd R the centroids of the plaquettes used in the correlator Eq. (3) 
lie at integer values of x\\. Figures 6 and 7 show the fields at the center of the Wilson 
loop as functions of R (the data for even R are averaged over plaquettes with centroids at 
x\\ = ±0.5). 

Near the quarks, the ratio of field strengths £ 1= i/£ J= i/2 is close to the Casimir ratio 8/3, 
for all R. In the region between the two quarks the ratio of fields falls below 8/3, reaching 
a minimum at x\\ = 0. The ratio is significantly reduced at large R. At R = 5, for example, 
£j=i/£j=i/2 = 1-97 ± 0.02 at T = 3, and 1.91 ± 0.08 at T = 4. Quantitatively similar results 
for the ratio of fields at a given R are found at all points in the plane perpendicular to the 
QjQj axis (see Figs. 8(a) and (b), described below). This suggests that there is a reduction 
in the color-electric flux in the region between the adjoint quarks. 

Some care is required in comparing the plaquette correlator data at fixed T for different 
i?, since the systematic errors in the T — > 00 extrapolation are expected to increase rapidly 
with R (the use of fuzzed Wilson loops to excite the QjQj pair becomes less efficient as R 
increases). The statistical errors for large R and T are too big to permit a reliable estimate of 
systematic errors in the extrapolation. Nevertheless, it appears that in some spatial regions 
the local adjoint field strength exhibits more pronounced screening than is indicated by the 
potential. This is presumably due to the fact that the potential averages the fields over all 
space. 

One naively expects that adjoint quarks become completely screened at sufficiently large 
i?, where it becomes energetically favorable for the flux-tube to fission, leading to the forma- 
tion of a pair of quark-gluon bound states ("gluelumps") [8,9]. One would therefore expect 
the fields in between the quarks to approach vacuum values at large R. A calculation on 
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larger lattices of the plaquette correlators at large R and T could shed light on this question. 
Unfortunately, the methods used here appear to be inadequate for this task. 

Results for Ej as a function of position x±_ perpendicular to the QjQj axis are shown 
in Figs. 8(a) and (b) (data was averaged over equidistant points along the two lines y and 
z from the center of the loop). The fields exhibit a "penetration depth" of a few lattice 
spacings. The fundamental and adjoint quark flux-tubes appear to have very similar cross- 
sections, despite appreciable screening of the fields in between the adjoint quarks at large 
R. 

This result has implications for the dual superconductor model of confinement [16], which 
is motivated by an analogy to compact Quantum Electrodynamics (QED), where confine- 
ment arises from vacuum condensation of magnetic monopoles [17]. In type II superconduc- 
tors a domain wall formed between normal and superconducting regions has negative surface 
energy, so that many normal regions are created, each carrying an elementary quantum of 
flux. In type I superconductors on the other hand a domain wall has positive surface energy, 
hence as few normal regions as possible are created [18]. If the dual superconductor picture 
is correct, then the fact that adjoint quarks form a single flux-tube, with a cross-sectional 
structure similar to the flux-tube formed by fundamental quarks, may suggest that the QCD 
vacuum is dual to a type I superconductor. Evidence for dual supercurrents associated with 
domain wall formation in SU(2) lattice gauge theory has recently been reported in Ref. [19], 
but the type of superconductivity could not be distinguished. 

These results also suggest a connection between the confinement mechanism in QCD in 
both three and four dimensions, and in three-dimensional QED. The fields generated by QQ 
pairs in several representations of three-dimensional SU(2) and U(l) lattice theories were 
calculated in Ref. [4]. The fields were found to be restricted to a flux-tube whose cross- 
section is approximately independent of the representation. It is interesting to note that 
Abelian magnetic monopoles are present in all of these theories. 
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FIGURES 



FIG. 1. Time-dependence of the estimated potential at R = 7 for J = 1 (•) and j = 1/2 (o). 
Statistical errors are shown when bigger than the plotted symbols. 

FIG. 2. Time-dependence of the ratio of fields £j = i/£j =1 / 2 a t the center of Wilson loops with 
R = 3 (□) and R = 5 (■). 

FIG. 3. Estimated potential versus quark separation R for J = 1 (•) and j = 1/2 (o). The 
j = 1/2 potential has been rescaled by C\j 2 = 8/3 (ci = 1). 

FIG. 4. Field strength versus position x\\ in the plane of the Wilson loop with R X T = 7 X 3, 
for j = 1 (•) and j = 1/2 (o). 

FIG. 5. Ratio of fields £j = i/£j =1 / 2 versus x\\ for Wilson loops with R = 3 (□) and R = 7 (■). 
In both cases T = 3. The dashed line shows the ratio of Casimirs. 

FIG. 6. Fields £j in the center of the Wilson loop versus R for j = 1 (•) and j = 1/2 (o). The 
data are at T = 3. 

FIG. 7. Ratio of fields £j = i/£j =1 / 2 m the center of the Wilson loop versus R. 

FIG. 8. Fields £j versus position xj_ perpendicular to the QjQj axis for j = 1 (A) and j = 1/2 
(o). The fields have been normalized to their values at xj_ = 0. The results are shown for (a) 
R = 3 and (b) R = 7. The data are at T = 3. 
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